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ABSTRACT 

In this paper we consider a massless scalar field, with a possible coupling ^ to the Ricci 
scalar in a D dimensional FLRW spacetime with a constant deceleration parameter g = e — 1, 
e = —H/H"^. Correlation functions for the Bunch-Davies vacuum of such a theory have 
long been known to be infrared divergent for a wide range of values of e. We resolve these 
divergences by explicitly matching the spacetime under consideration to a spacetime without 
infrared divergencies. Such a procedure ensures that all correlation functions with respect to 
the vacuum in the spacetime of interest are infrared finite. In this newly defined vacuum we 
construct the coincidence limit of the propagator and as an example calculate the expectation 
value of the stress energy tensor. We find that this approach gives both in the ultraviolet and 
in the infrared satisfactory results. Moreover, we find that, unless the effective mass due to 
the coupling to the Ricci scalar ^-R is negative, quantum contributions to the energy density 
always dilute away faster, or just as fast, as the background energy density. Therefore, 
quantum backreaction is insignificant at the one loop order, unless ^-R is negative. Finally 
we compare this approach with known results where the infrared is regulated by placing the 
Universe in a finite box. In an accelerating universe, the results are qualitatively the same, 
provided one identifies the size of the Universe with the physical Hubble radius at the time 
of the matching. In a decelerating universe however, the two schemes give different late time 
behavior for the quantum stress energy tensor. This happens because in this case the length 
scale at which one regulates the infrared becomes sub-Hubble at late times. 
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1 Introduction 



It is by now well established that on the largest scales, the universe is well described by a 
homogeneous, isotropic and spatially flat metric, given by [1] 

g^ydx^dx" = —dt^ + ait^dx ■ dx. (1) 

Here a is the scale factor, and the derivatives on a deflne the Hubble parameter, H and the 
deceleration parameter q (in this paper we shall use the equivalent parameter e) 

H^- ; g = -l-^ = -l + e. (2) 

If e > 1 the expansion of the universe is decelerating, if e < 1 it is accelerating and if e = 1 
it is curvature dominated. In a typical cosmological model, we couple this geometry to a 
certain perfect fluid, with energy density = pi,{t) and pressure pf, = Phif)- The dynamics 
if the scale factor are then governed by the Friedmann equations 

^ - (3) 
H = -4:TcGN{pb + Pb) , 

where Gn indicates Newton's constant and for simplicity we assumed flat spatial sections. 
This is consistent with the current cosmological measurements [2j, which give a very large 
lower limit on the radius of curvature Rc of the Universe, Rc > 22h~^ Gpc {Rc > 33h~^ Gpc) 
for a universe with positively (negatively) curved spatial sections, where h = 0.705 ± 0.013 
is the rescaled Hubble parameter (in units of 100 km/s/Mpc). If we now assume that the 
cosmological fluid obeys an equation of state 

Pb = uibpb, (4) 

with w constant we immediately flnd that in such a case 

3 

e = 2 ~'~ ~ constant. (5) 

The equation of state (jll) is in practice obeyed for many types of fluids. For example radiation 
corresponds to Wb = 1/3 and e = 2, non-relativistic matter has Wb = and e = 3/2, the 
cosmological constant has Wb = —1 and e = 0. If e is a constant, we can flnd the Hubble 
parameter and the scale factor function of time 

if(t) = ^^- ; a{t) = {l + eH,t)\ (6) 
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with Hf) a constant. 

The behavior of a massless scalar field with a possible coupling to the Ricci scalar on such a 
background geometry has been extensively studied in the literature [3] [1] [5] [6j . An interest- 
ing observation is that the kinetic operator of, for example, the graviton can be related to 
the kinetic operator of such a massless scalar field. Therefore the propagator of the graviton 
can, apart from some tensorial structure, be written in terms of propagators of massless 
scalar fields [8J [9J . Similar observations can be made for vector fields [20] or anti-symmetric 
tensor fields [31j. Therefore an understanding of the physics of the massless scalar field can 
be translated to those fields as well. 

In particular many studies concern in particular the special limit of e (de Sitter) [7j. 
Because of the nature of the background space-time, the scalar field has non-trivial quantum 
properties. In particular, due to the presence of a horizon, there is particle production. The 
particles are created mostly in the infrared, and it is this effect that leads for example to 
the creation of primordial density fluctuations. Also it is the basis for many works of the 
quantum backreaction on the background space-time [10] [H] [12] [13] [E] [I5] [16] . However, 
it was long ago realized that when one chooses the vacuum to contain only purely positive 
frequency modes (the so called Bunch-Davies vacuum), the expectation of the two point 
correlator for this vacuum diverges in the infrared, for a large - and physically relevant - 
range of e [5]. In particular, if the massless scalar fleld is minimally coupled, there are in- 
frared divergencies for all e < 3/2. The presence of these divergencies implies that in the 
infrared, the Bunch-Davies vacuum does not describe a physically sensible state. One possi- 
ble solution to this problem is to assume that the spacetime manifold is spatially compact, 
for example a torus T^~^. This approach effectively introduces an infrared cut-off, when 
the sum over the modes is approximated by an integral [VT}. In Ref. [18] the propagator is 
explicitly constructed using this regularization for any constant e space. 
A different approach is to choose the mode functions such that the super-horizon modes are 
less singular than they would have been in the Bunch-Davies vacuum [SI S]. Because only 
the super-horizon modes change there would be no effect on the Hadamard short distance 
behavior of the propagator. The time dependence of the mode functions is determined by the 
scalar fleld equation but their initial values and the initial values of their flrst time derivatives 
can be freely specifled. Thus we can choose the initial values for the super-Hubble infrared 
modes to be in some infrared flnite state. Such a choice would ensure then that there would 
be no infrared divergence, either initially or at any later time |19j . 

In this paper however, we propose a third regularization of the infrared. Instead of splitting 
the ultraviolet and the infrared sector, one could also consider the matching between the 
Bunch-Davies vacuum in an infrared safe space-time and the space-time one wishes to study 
for all modes. Also in this case the initial state does not lead to infrared divergences and thus 
also the flnal state will be safe. In this approach the ultraviolet mode functions will differ 
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from the ultraviolet Bunch-Davies mode functions, but the ultraviolet divergent structure 
will not change. The advantage of this approach over the previous two is that one can quite 
easily envision something like this to be realized in nature. For example if inflation was 
preceded by a radiation dominated epoch, we essentially have the scenario described above, 
if the field was in its Bunch-Davies vacuum during the radiation epoch and if the transition 
between the two geometries is fast enough. A disadvantage is that the sudden matching 
between the two geometries causes the Ricci scalar to change discontinuously. This causes a 
burst of particle creation, but we shall see that this does not significantly influence our main 
results. 

There are other regularization schemes proposed in literature. Here we mention a recent 
proposal by Parker [20], where it was argued that the infrared should be regulated by sub- 
tracting the adiabatic vacuum contribution. The proposal has been worked out to some 
detail in Ref. [21]. We feel that that this scheme is not well motivated, since the infrared 
sector of constant e spaces strongly breaks adiabaticity, which is the assumption made when 
constructing adiabatic vacua. Yet another possibility is to subtract the vacuum of comov- 
ing observers [22]. This procedure works well for de Sitter space p2], but it remains to 
be seen whether such a procedure can be successfully implemented in more general FLRW 
spacetimes. 

In this paper we want to study the role of the quantum infrared fluctuations generated 
by the expansion of the Universe for a massless scalar field. In particular we would like to 
understand whether the energy density due to these fluctuations could ever become relevant 
for the evolution of the universe. For a proper understanding of such a question, it is im- 
portant to compare different regularization schemes for the infrared. In particular we expect 
that in an accelerating space-time, results will not depend that much on the regularization 
chosen, since the details of any such scheme will grow quickly to super-Hubble scales. The 
reason is that in accelerating space-times, physical length scales grow faster then the Hubble 
radius. In a decelerating space-time however, the opposite is true. Specific cases need to be 
studied to understand this issue in detail. 

In section [2] we discuss the construction of the scalar field mode functions and its infrared 
properties. In section [3] we construct the new vacuum by matching it to an infrared finite 
Bunch-Davies vacuum. The associated propagator is constructed in section HI We discuss 
some properties of the stress energy tensor in section O and we calculate the ultraviolet 
contribution to the stress energy tensor in section and consider its renormalization. In 
section [7] we consider the equation of state of the quantum fluid and study the scaling of 
the quantum energy density. In section [8] we compare our results with the results obtained 
in Ref. [18] using a cut-off procedure and finally in section [9] we summarize our results and 
conclude. 
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2 A Massless scalar field 



We wish to study the behavior of a massless scalar field in the geometry ([T]). The action 
in D dimensions is given by 



(7) 



where ^ is a constant coupling of the scalar field to the Ricci scalar R. The equations of 
motion for (f) are 

(□ - iR)<P = 0, (8) 
where in our background metric the scalar d'Alembertian reads, 



92 1 52 d 



Now we use (spatially) Fourier transformed field 



k) = / ci^-^xe~*'=-Xa;) 



to obtain the equations of motion for 



+ {D- l)H<p +—(P + iR<p = Q {k 



(9) 



(10) 



(11) 



where = d(j)/dt. We expand our solution in (time independent) creation and annihilation 
operators 6^ and h 

;f)(t,k) = ^{t, k)h{k) + ^*(t, k)h\-k) , (12) 

where the spatial homogeneity of the mode functions ip{t, k) implies that they are a function 
of the modulus of fc, and h{k ) is the annihilation operator that annihilates the vacuum 
h{k)\VL) = 0. The normalization is fixed by requiring the canonical commutation relations 



h(k),h\k') ={2Tx)°-H°-\k-k') 
'4>{t,k), a{tf''4>^{t, k' )j = «(27r)^-i5^-i(fc -k'), 

which implies that the Wronskian of the mode functions ip is given by 

ipit, k)tp*{t, k) - tlj*{t, k)^{t, k) = ta^-^ . 



(13) 



(14) 
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The mode functions ip obey the equation 

where we used that in the constant e geometry under consideration [23] 

R= {D -1){D -2e)H'^. (16) 
We write z = jiz^jjj^- If we assume that e is constant as in ([5]), we obtain 

where 

^ / D- 1-6 x2 (D-l)(D-2e) ^ 
V 2(1 -e) ; (1 -e)' 

1 {D-2e)f D-2 



4 (1-el 



■(P-I)e 



such that V = 1/2 for a conformally coupled scalar, for which ^ = = (-D — 2)/[4(D — 1)], or 
when the universe is radiation dominated, such that e = D/2. Equation (ITTj) can be easily 
related to Bessel's equation and we find the solution to be 

i){t, k) = a^~^ (a{k)u{t, k) + l3{k)u*{t, k) 

(19) 



u{t^ k) 



m-e)Ha \(\ 



where h'^^ denotes the first Hankel function. The normalization is chosen such that the 
Wronskian condition (fT^ implies that 

= 1. (20) 
The time-ordered (Feynman) propagator for a given vacuum state is now given by 



i^{x;x') = {n\e{t - t')(i){x)(f){x') + e{t' - t)(i){x')(i){x)\n) 

(27r)^-i 



fjD-^U ^ ^ / N f21) 

^H^-')(e{t-t')^{t)^*{t') + e{t' -t)^{t')^*{ty ^ ' 
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The vacuum \Q) is defined by b{k = 0. This vacuum state is clearly not unique, since 
the commorients a and f3 are not uniquely determined by fl20l) . The particular choice 



a = l ; (3 = 0. (22) 

corresponds to a field in a state \Q) known as the Bunch-Davies vacuum [3|. This choice 
minimizes the energy in the ultraviolet modes and is therefore often considered as a reason- 
able choice for the vacuum. By considering the small argument limit of the Bessel function 
(see Eq. (8.402) of Ref. [23|) 



ra=0 



and using that (see Eqs. (8.403) and (8.405.1) of Ref. [25) 



we find that the small k contribution to the integrand of fl2T|) in the case of the Bunch-Davies 
vacuum behaves as 

iA(x;x)BDOC J k^-^-^''{l + 0{P))dk, (25) 



where we took z/ > in Eq. ( 1181) . We thus see that the propagator of the Bunch-Davies 
vacuum is infrared divergent for all [5], |18] 



D-1 

V > ^— . (26) 

For example if = and D = A (we do not need dimensional regularization in the infrared) 
this implies that there is an infrared divergence for all e < 3/2. From ([5]) we find that this 
implies that the pressure of the fiuid driving the Universe's expansion is negative. In the 
more general case, where ^ 7^ 0, we refer to figure ([1]), where we plotted (in D=4) the regions 
in the e — ^ plane where the requirement fl26|) is met. We can qualitatively understand this 
picture, by realizing that the coupling ^ acts as an effective, time- dependent mass term for 
the scalar field proportional to C,R = 6.^(2 — e)if^. Thus for e > 2. Thus we have an infrared 
divergence if this mass term is 'sufficiently' less then 1/6 (where the precise meaning of 
'sufficiently' follows from (!26|) . 
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Figure 1: The shaded regions indicate the regions where > |, with v given in f|T8l) and 
we choose D = 4. When this requirement is met, the coincident propagator f l25|) is infrared 
divergent. The regions are bounded by the curve ^ = and the dotted asymptotes are 

e = 2, ^ = 1/6 and i = 3/2. 
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3 Matching 



To resolve the infrared divergences, we shall consider a geometry where v = 1/2 at times 
t < t, while z/ is an arbitrary half integer for t > t. Notice that u = 1/2 corresponds to a 
radiation dominated universe (e = D/2) if ^ = 0, or a universe whose expansion is driven by 
a conformally coupled scalar = 4^1^ )• The restriction that u after the matching should 
be half integer comes from technical considerations. The calculation simplifies enormously 
in this case, since the mode functions can now be written as finite sums. Moreover we shall 
find that in the end of the calculation we can resum these finite sums and analytically extent 
the result, obtaining thus an answer, which then will be valid for arbitrary u. 
We match these two geometries, such that the scale factor a{t) and the Hubble parameter 
H{t) are continuous at the matching. This requires that also the mode functions and their 
first derivative need to be continuous at the matching [25j|. It is clear from fl26l) that before 
the matching, in the u = 1/2 spacetime, we can freely choose the Bunch-Davies vacuum. 
Since there is no infrared divergence initially, it will not be there at any later time. 
Quantities before the matching we indicate with a subscript 0, while quantities after the 
matching do not have a subscript. Thus for t < i we have eo and for t > i we have e. The 
mode functions f|T9|) are written in terms of the following solutions 



1 ik 

uo{t < i, k) = =e(i-^o)ff- 



4:{l-e)Ha " \{l-e)Ha) (27) 



1 „ 1 v4 r(z/+i + n) / -2ik 



^n!r(z/+i-n)V(l-e)^ra 

where the last line follows from the fact that i/ is a half integer [24j. We choose the initial 
vacuum to be Bunch-Davies, and therefore the mode functions are 



ipoit <i,k) = ^Uo(t,k) 
M > i, k) = a^-^ (a{k)u{t, k) + (3{k)u*\t, k) 



(28) 

We shall now fix the coefficients a and (3 by requiring that at the matching a and H are 

dt 

Mi k) = tp{i, k) 



continuous and thus ip and 4'?/' should be continuous 



d 



dt 



d 



-Mt,k) ^=-M,k) 



t=t 



dt 



(29) 

t=i 
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These conditions imply that 



a = ial u*Uq — u*Uq 



) 



t=t 




— I — 



{l-e)np\ V{u+\ + n) / 2ik 



2k ) n\r{u+^-n)\{l-e)p 



) 



— n 



(30) 



^-if+a^e P (i-€)(i-€o) 




{1 — e)np\ T{iy+^ + n) / — 2ifc 



2k )n\r{u+\-n)\{l-e)p 



) 



— 71 



n=0 



where we defined p = a{i)H{i). Notice that the series expansion used for the Hankel functions 
is also valid also for negative values of the argument, so we do not run into any problems 
due to the branch cut of the Hankel function for negative real arguments. 

4 The Coincidence Propagator 

Now we shall calculate the propagator, using the state \VL) defined by fl30l) . In particular, we 
shall only calculate the propagator at coincidence. As we shall see this is sufficient for our 
present purpose, calculating the one loop corrected stress-energy tensor. The coincidence 
propagator is given by 





(31) 
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iA(x: x) 



dkk 



D-3 



Using the expressions from the previous section we obtain 

(4vr)Vr(^ 

r(z/ + i + n)T{u + i + m)T{u + i + q)T{u + i + r) 



,r,m,n=0 \ ^ ^ ' ^ ^ ' ^ z ^ z ' 



1 + -(1 - e)(g - n)f + -(1 - eYnq^- (-1) 



A;2 



(32) 



+ e(i-^)-f^-^ p ^ - -(1 - e)g^ - -(1 - eYqn 



+ e( 



4 



;i - e)p 



2ik 



[1 - e)aH 



n\m\q\r\ 



One important property of this expression is that the integral converges in the infrared, or 
— >■ 0. Of course this was expected by construction, but it can be shown exphcitly. For 
example the quadruple sum in (1521) evaluates in some explicit cases to 



E 

q,r,m,n=0 



E 

q,r,m,n=0 

-i 

E 

q,r,m,n=0 



(33) 



=7/2 \ p / \7 1 aH' 



4.1 The Ultraviolet 

As long as aH ^ p, all UV divergencies can be easily seen to come only from the first line 
in square brackets in Eq. fl32|) . In fact the second line has three types of contributions: 
UV divergent terms, IR divergent terms and a term that is logarithmically divergent in 
both the UV and the IR. Now we know that the IR divergent terms will cancel against 
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the IR divergent terms coming from the second and third hnes in the square brackets. So 
we shall simply drop all IR divergent terms that we encounter (except for the logarithmic 
divergence). By explicitly calculating the first 3 terms of the sums, we find that the UV 
divergent contributions are 

^A{x■, x)uv = f'l^ I (k^-' - ^^k^-%1 - era - Au^))dk (34) 
(47r)~r(-^) J ^ o J 

The first (quadratically divergent) term is a scaleless integral, which we can automatically 
subtract in dimensional regularization [26] [27]. The second contribution is logarithmically 
divergent. Since we know that the infrared divergent part should drop out in the final answer, 
we write the integral as 




(35) 



and we drop the first, IR divergent, integral. Now we obtain 



. ^) [_ _ I ^ (^) + 0(Z, - 4)] . (36, 

where one should note that the prefactor (1 — e)^(z^^ — |) is still D-dependent, as can be seen 
from f|T8l) . and we introduced a renormalization scale yU. The k^ dependence of fl36l) should - 
and does - drop out when we add the IR contribution. Notice that A(a;; x)\jy is exactly the 
same as the (divergent) UV contribution one would have obtained using the Bunch-Davies 
vacuum [18J. This was to be expected, since the ultraviolet of the theory is only sensitive to 
local physics and therefore does not 'see' the matching. 



4.2 The Infrared 

We now focus on the second line in square brackets of fl32|) (the third line is simply the 
complex conjugate of the second line). As long as aH ^ p, the contributions are UV finite, 
so we can put D = 4. We shall comment below on the special case aH = p. The integrals 
are of the form 

j k-'^e^'^dk, 

where n is a positive integer and a < 0. We split the integral again in two ranges as in 
( !35|l and take the limit /cq ^ 0. The lower integral will be divergent, but these divergences 
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exactly cancel against the infrared divergences coming from the second line of fl32l) . The 
upper integral is given by {'^[a] > 0) 



k-''e'°''dk = kl-^'Eni-iako) 



(37) 



ko 



where n is an integer and 



dt 

tp 



denotes the exponential integral. For integer index, the exponential integral has the following 
expansion around zero 



\n— 1 



Eniz) 



Tin) 



ipin) — ln{z) 



(k-n + l)k\ 



(38) 



where ip{z) = (d/dz) ln[r(z)] denotes the digamma function. Since we know that all negative 
powers of ko cancel, we obtain the following leading order result, which is obtained by 
expanding (j37j) around ko = 0, 

(±ia)"-i 



ko 



Tin) 



[ln{Tiako) - ^{n)] + (IR - div) + O{ko) 



(39) 



where ip{z) = {d/dz)ln[T{z)] denotes the digamma function. 

We use this expression to evaluate the third and fourth line in fl5^ and obtain 



iA{x;x)iR 



r(z/ + i + n)T{u + i + m)T{iy + i + q)T{u + i + r) 



47r2 ^ r(z/ + i - n)T{iy + i - m)T{iy + i - q)T{iy + ^ - r) 

n,n,q,r=0 ^ ^ ^ ' ^ z ^ A ' 

ln(|a|A;o) ~ 'i^ija + ?T, + g + r + l) ln(|a|A;o) — '^{m + n + g + r) 



X 



r^(i-C) 

T37(-(i-0) (7-1 



T{m + n + g + r + r 

Q C"^ /' \m+r \ 



T{m + n + q + r) 



nlmlqlrl 



(40) 



where a 



;i-C"^) andC = p/{aH). 



{l-e)aH ' 

Quite remarkably, the sums in Eq. (HOl) can be (almost completely) performed, resulting 



m: 



iA{x;x)iR 



H\l-ef 
87r2 



n=l 



n=2 



(41) 
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where Cy is a slowly varying function of Its precise value is not important, as it can be 
absorbed in the by a finite counterterm. 



4.3 The Full Propagator at coincidence 

The full propagator at coincidence is given by the sum of the infrared and ultraviolet con- 
tributions. Adding the infrared fBTl) and ultraviolet fl5B]) contributions we obtain, 



4 

iy-3/2 



/i 



D-4 



D-4: 

2^ ^ 2n ^ 2n 



^ + iln 
2 2 



^+1/2 ^2n' 
2 

n=l n=l 



(1 

lla) ■ 



1x2 



(42) 



We see indeed that, as promised, the two logarithms of ko coming from the infrared and the 
ultraviolet have canceled. From the result ( l42i) we see immediately that when |(^| <^ 1 the 
leading order contribution to iA{x;x) is of the form, cx H'^(^~'^'^, while when \(\ ^ 1 the 
leading order contribution goes as oc H'^C,'^^^^. 

The two series in (1421) can be resummed, resulting in: 



iA(x; x) 



167r2 



2;U 



D-4 



D-4 



- Cj, + In 



2^2 



1-- 



c 



+2C-C' 



^1-2U 111 

X 2FAI, u--; Z/+-; — 



z^-t V 2 2'C 

C'+'" ^ 3 5 



ln(l-C^)-^x,F,(l,.+|;.+ ^;C^ 



(43) 



We would like to interpret this result as the analytic extension of Eq. fH2l) to arbitrary 
(complex) u. In order to do this uniquely we need to specify the Riemann sheet of both the 



^We were unable to evaluate the constant c^. The first few values are given by quite simple expressions. 
For ly = {3/2, 5/2, 7/2, 9/2, 11/2} we have = ij{2iy - 2) + {1/2, 1/3, 1/2, 2/3, 137/168}. 

^Mathematica represents the answer in terms of the Lcrch transedent, LerchPhi[z, s, a], which is the 
following generalization of the Riemann C. function, LerchPhi[z, s, a] = X^J^o z"' /{i^ + ctY ■ For our purpose 
it is convenient to express the Lerch transedent in terms of the Gauss' hypergeometric function. 



LercliPhi[z, 1, a] 



2^1(1,0;; a + l;z) 
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logarithm and of the hypergeometric function in the second hne for |C| < 1 and in the third 
hne for |C| > 1 (recall that the logarithm has a branch cut along the negative argument and 
that the hypergeometric function has a branch cut running along positive arguments z > 1). 
It turns out that the following cut prescription uniquely specifies the analytic extension to 
arbitrary u: 



ra=l 



n 



u+1/2 



E — 



n=l 



ICI<i 



ICI>i 



1 



1 



-In 



1 



/■l-2i/ /II 

x.F.(l,.--;. + -: 



1 



2' C±ie 



C / 3 5 

J— X 2Fii^l,--iy;--iy;C) +7rtan(7rz/) 

3+2u 



(44) 



-In 



1 - 

1/e 



2^ ^+2 

^2.+l ^ / 1 1 1 



+ 7rtan(7rz/) , (45) 



where e > is an infinitesimal parameter. Notice that the results in ( |44] - I45l) are independent 
on the cut prescription. We are now ready to write down a complete expression for the 
propagator at coincidence. In the accelerating case (e < 1, C = p/{ciH) < 1) from Eqs. (H3l) 
and (HH) we thus have, 



iA{x; x) 



HHl-eY / 2 1 



167r2 



2/i 



D-4 



Cu + 7rtan(7rz/) + In 



C^^" „ 3 , 5 



+ 2C-C^-El^^x^^i(l'i±-5^±-'C^ 



(46) 



Notice the symmetry v —>■ —v in all C dependent terms (although, one should keep in mind 
that we assumed > 0). This symmetry has been observed when the infrared is regulated 
by placing the Universe in a comoving box [18]. On the other hand, in the decelerating case 
(e > 1, C > 1) we get from Eqs. (|13]) and fjltj) . 



r2n .^2 ^ 



167r2 



2/i 



D-4 



L)-4 



— Cj, + TT tan(7rz/) + In 



47r/i2 



l-e)2/72(i + ^)2 
+ 2C-C^-E\^X2F,(l,--±.;-±.;-) 



(47) 
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Just like the accelerating case (H^ . this expression also exhibits a z/ — —v symmetry in all C 
dependent terms. It is worth noting that in the limit when i/ is a half integer, z/ ^ + 1/2 
[N = 0, ±1, ±2, ..), both coincident propagators fH6l) and P7|) are finite. In the special case 
= 0, both are simply zero, while for all other cases the simple pole in the tangent cancels 
against the simple pole in the hypergeometric function. 

From Eqs. (HB1 - [17|) we can easily extract the leading order late time (a oo) behavior 
of the coincident propagator in the accelerating case {aH oo, such that ( ^ 0) and 
decelerating case {aH 0, such that ( ^ oo): 



In the special case when u = 3/2 the dependence on ( in the accelerating case is logarithmic. 
4.4 Numerical check of the analytic extension 

To test whether the procedure described above to extend the propagator at half integer u 
to all values of u, we have performed some numerical checks. To get a numerical result 
for the propagator, we use the mode functions in terms of the Hankel function, as in flTI?]) . 
From those we can easily get the general expression for a and (3 from the first lines of (130|) . 
When evaluated numerically Eq. ( 12T|) then gives an integral which, after the appropriate 
regularization, corresponds to (HGl) . This is, of course, true up to the 1/{D — 4) term, 
which is removed by renormalization. There are two ultraviolet divergences: a quadratic 
one and a logarithmic one. We can easily read off both divergences from (IMI) . We know 
those expressions for the divergences are also valid for non half integer u. The quadratic 
divergence we can simply subtract from our numerical result. In order to correctly regulate 
the logarithmic divergence, we choose an ultraviolet cutoff A for our numerical integrals and 
realize from flM|) that we must identify 



To smoothen the numerical procedure, we also choose an infrared cutoff kg. Apart from the 
1/(D — 4) term, the numerically integrated result should now correspond to (H6l) for all values 
of z/, up to the undetermined constant Cj, in the analytic result P6l) . We fix this constant 
for a certain fixed, but arbitrary set of parameters. As an example we show the result for 
u = 1.7 in figure [21 The main plot shows the behavior of the coincident propagator divided 
by the Hubble scale A{x; x)/ H"^ versus conformal time defined by 





ln(A) = ln(a/iv/vr) + 1 - -7^;. 



(49) 



1 



(50) 



{l-e)Ha 
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For the matching point we choose p = 1. Since e = 0.1, we thus find that fj ^ —1.1. The 
inset in the figure shows the difference between the numerically calculated propagator and 
the analytic solution fH6|) . We thus find that, up to a numerical noise of the order of 1% 
or less, the two results agree. Such a noise level is not unreasonable given the complexity 
of the numerical integrations. We have also checked several other values of z/, which give 
similar results. Albeit we do not have a rigorous mathematical proof, these numerical results 
strongly suggest that the analytic extension to all u made in fHB]) and flTTI) is correct. 



A(x;x)/H- 




-1.0 -0.8 -0.6 -0.4 -0.2 0.0 



Figure 2: The main plot shows the coincident propagator ( H6l) . rescaled by versus confor- 
mal time ( l50l) . The inset shows the difference between the analyticially calculated result and 
a direct numerical integration of (!2T|) . using the mode functions and the a and f3 coefficients 
in terms of Hankel functions (fT9|) . The parameters chosen are p = 1, a = 4, u = 1.7, e = 0.1, 
A = 10^, ko = IQ-^. We find for this case that ^ tlj{2u - 2) + 0.348. 
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5 Energy momentum tensor 

The energy momentum tensor is defined by 



By varying the action ([7]) we obtain 

T^u = - + i{R^, - ^g,.R)^' - e(V^V, - g,.n)(p\ (52) 

We take the trace to obtain 

T% = ~^^[{d'^<l))d,<P + + iD- l)^n<p' (53) 

and, using the equation of motion ([8]) for 0, we find that the one loop expectation value is 
given by 

T, = {Q\T%\Q) = ((D - 2) - A{D - l)^)aA{x; x) . (54) 

Since any quantum correction will respect the symmetries of the underlying spacetime, we 
know that we should be able to write (we use a subscript q to indicate we are considering 
the quantum corrections to the stress energy tensor) 

D-l 

and then we can calculate pq and Pq from Tq using the covariant stress energy conservation, 
which we write as ^ 

— Pq + Dpq = Pq-{D - l)pq = -Tq . (56) 

Now an important question is whether the energy density in the quantum corrections can 
dominate over the energy density in the classical background. The background Friedmann 
equations ([3]) tell us the scaling of energy density 

where Wb denotes the equation of state parameter of the background and pb = Pbit). Since 
the background energy density is responsible for the expansion of the Universe, from ([5]) we 
have that 

Wb = — = -l + -e. (58) 
Pb 3 
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Furthermore, we know from the conservation equation fl56l) . which of course holds for any 
stress energy tensor, that the energy density of the quantum contribution scales in general 
with the scale factor as 

where Wq = Pq/ Pq is the the equation of state parameter for quantum matter contribution 
and Pq is the energy density at the matching, which is typically small when compared to the 
background density at the matching, pg ~ ~ {H/Mp)'^pb. Here Mp = (SttGat)"^/^ ~ 
2.4 X 10^^ GeV is the reduced Planck mass, H = H{t) is the Hubble parameter at the 
matching. For simplicity we took the scale factor at the matching to be equal to unity, 
a = a{t ) = 1. 

From Eqs. (15714591) we see that 

Pi = Pla^i-'t-^,) ^ f JLVaH-^^--'.) , (60) 



Pb Pb vMj 



which implies that the quantum contribution to the stress energy grows with respect to the 
background contribution whenever 

Wq<Wb. (61) 

Hence the principal task of our investigation is to find out whether this condition is ever 
met. If the answer is affirmative then - in the light of Eq. fl60|) - we would be lead to 
the conclusion that (after a sufficiently long time) the energy of quantum fluctuations would 
eventually dominate the energy of the background, leading to a strong quantum backreaction. 
Such a strong backreaction might then signiflcantly change the evolution of the Universe. Of 
course this would then in principle ruin our Ansatz that e is constant, and thus one should 
be careful in interpreting the result. If the backreaction is large, the critical time tc after 
which quantum contributions start to dominate can be estimated to be, 

tc = t\ — \ ~H~^] {Wq<Wb). (62) 

^Pq' ^ H ' 

If however Wq > Wb, the energy density due to the quantum contributions will be negligible 
at all times. But before we address this question, we must renormalize the stress-energy, 
which is what we do next. 
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6 Renormalisation 

The only divergence arising in the one-loop stress energy tensor is induced by the divergent 
part of the coincident propagator (143!) , which we write in the form 



where we took account of Eq. ( 1181) . Now from the trace equation (IMI) we get for the divergent 
contribution to the stress energy trace, 

[{D-2)-4{D-imD-2e){D-l-3e)e^,^-'H' 



?7r2 D-A 167r2 



(19-23e+6e^) - 6^(15-17e+4e^ 



where we took account of UH'^ = -{dt + (£> - l)H)dtH'^ = 2t{D - 1 - 3e)H^. 

It is known that this theory can be renormalized using only one counterterm, proportional 
to E?. Indeed, taking a functional derivative with respect to g^'^ of the B? counterterm action 
results in [181 
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(ct)^, = --^^ j d^'x^aE? = a(4V^V,i? - Ag^^UR + g^,R^ - ARR^,) . (65) 

Multiplying by g^'^ results in the trace 

(ct)^'^ = a(^-A{D-l)nR+{D-A)R'^^ 

= a{D-lf{D-2e)H^i^ - 8e{D-l - 3e) + (D-4)(D-2e)) 

= -432a(2-e)(l-e)ei/^ + 36a(4-34e+35e2-8e^)(D-4)i7^ 

+ 0{{D-Af). (66) 

Upon comparing this with Eq. (!64l) we can read off a that renormalises the theory, 

" = Tl52;^D^ + "^' 

where is an arbitrary but finite constant. One can easily show that this choice of a renders 
all components of {Ti^)q finite. Combining the two contributions (IM|) and (|66|) . we get the 
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following finite expression coming from the UV divergent contribution, which we indicate 
with a superscript (1) 

T,<" = (T,)di. + (ct); 

= ^i^5^£?^^((l-60 + 4^)- 432«,(2-e)(l-.).//V (68) 

In the following section we shall consider the contribution to the trace of the stress 
energy tensor coming from fH3l) to study the late time behavior of the quantum fluid. This 
contribution we shall indicate with a superscript (2). 



7 The equation of state of the quantum fluid 

In this section we calculate both the pressure and energy density of the quantum fluid 
that results from one loop scalar fluctuations in our model. The general procedure is quite 
straightforward (albeit somewhat tedious): we first need to act with the d'Alembertian on 
the ultraviolet parts of the coincident propagator given in Eqs. ( 1461 - I471) to get the trace of the 
stress energy tensor Tq as defined in Eq. (|54l) . To get the full Tg one has to add the finite 
contributrion Tg^^ in (1681) that remained after renormalisation. When we have obtained the 
full stress energy tensor, we can solve the conservation equation fl56l) to obtain the quantum 
energy density pg and quantum pressure Pg, from which we obtain the quantum equation of 
state parameter Wg = Pg/pq- To facilitate this procedure we give a list of useful formulae in 
Appendices A and B. 

To illustrate the procedure, we shall first calculate the contribution to pg and Pq coming 
from the finite remainder from ultraviolet fluctuations (l68l) . We denote these contributions 
by p^g^ and p^g\ 

To calculate p^g^ we need to integrate Eq. fl56l) . which we can do with the help of Eq. (11351) 
in Appendix B. The result is, 

p« = -^^^|^^T^((l-60 + 46) + 108a,{2-e)eH' (69) 

The contribution to the pressure is then simply. 



= P^^{-^ + t') (70) 



Ti'^'+P? f (l-60(2-e)2i74, , ..^,.^,.no ^rA^ .4 
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Of course, the one loop infrared contribution needs to be taken into account for a complete 
answer. We shall now consider the general infrared contributions to the quantum energy 
density and pressure. We first consider the accelerated and then the decelerated case. 



7.1 Matching onto acceleration (e < 1) for a general u 



In this section we consider in some detail the general accelerating case (e < 1, C < !)• The 
decelerating case we shall consider in section 17.21 The contribution to the trace of the stress- 
energy tensor coming from the UV divergent terms has been calculated in fl68|) . Our starting 
point for the remaining terms is the coincident propagator fH6|) . Now making use of Eq. (11301) 
in Appendix A the contribution to the trace of the quantum stress-energy tensor flM|) from 
the remaining terms becomes: 



T(2) 
1 



:i-60^(l-e)(2-e)g^ 
167r2 




2 V(l-e)2i72(i + ^)2y 2 



3±2z/ 



X 2F1 (l,-±v--± V- j + - tan(7rz/) 



3-5e 2(2-3e 

— + 



1-e 



1+C (1+0- 



+ 2(l+e)C-C' 



(71) 



5^[(3-56)-(l-e)(3±2z.) 
± 



c 



3±2!/ 



+ 2(l-e- 



-c 



2^2 



Recall that by assumption u > 0. The u = 1/2 pole of tan(7rz/) is harmless since in this case 
^ = 1/6, such that the whole expression vanishes as the result of a vanishing prefactor. In 
Eq. ( I7T1) we made use the derivative of the hypergeometric function. 



C—^y<2Fiil,a;a+l;e) 
dQ 2a 



c 



2a 



The total stress energy trace is then obtained by adding Eqs. (!68l) and (ITTIl . 
A careful look at Eqs. (I7T1) reveals a quadratic divergence in Tg at the matching time t = 
t, at which ( ^ 1. This divergence is also present in the energy density and pressure 
calculated below in Eqs. fl75]476|) . and shows up as the logarithmic divergence in the coincident 
propagator fj43|) . This divergence occurs because of the sudden matching procedure assumed 
in this work, where e jumps suddenly from e = 2 to an arbitrary e, which also implies the 
sudden change in the Ricci scalar and the scalar field mass parameter. Of course, in any 
physically realistic situation any of those parameters will change smoothly. These sudden 
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changes induce a change in the quantum contribution by an infinite amount at the matching. 
An attempt to renormahse it away by including it into the 1/{D — 4) subtraction in Eq. fl34l) 
would result in divergent contributions away from the matching, which would make matters 
only worse. But let us try to understand how bad the matching divergence actually is. 
From Eq. fl7T]) we see that its contribution to Tq is of the order ~ H^/[^ — aH / [aH)]'^ . By 
taking t — i = 5t small, and requring that pg < = 3MpH'^, where Mp = (SnG)^^^"^ is the 
reduced Planck mass, we get that pg < ph implies 6t > — e\Mp], which is of the order 
the Planck time. Based on this estimate we expect that the matching divergence will be 
regulated whenever a sudden matching is replaced by a smooth matching, in which e at the 
matching changes over a period of time that is much longer than the Planck time. 
Before we calculate the quantum energy density, we take one more look at the form of 
Eq. fl71l) . All terms containing u dependent powers of ( in fl71l) can be written as the 
following sum. 



:i-60^(l-e)(2-e)g^ 
167r2 



EE 

± n=0 



6e 



3±2u+2n 



+ (3-5e) - (l-e)(3±2z/+2n) 



3±2u+2n 



(72) 

The contribution of this sum to the energy density can be, up to a z/ dependent integration 
constant cte(z/)(^^if^, determined by integrating Eq. fl56|) . which can be performed by making 
use of Eq. fll39p . The result is 



:i-6o^(i-6)(2-6)g^ 

167r2 



EE 

± n=0 



6e 



3±2z/+2n 

^3±2i/+2n 



+ (3-5e) - {l-e){3±2iy + 2n) 



'l-Qif{2-e)H^ 



(l-e)(-l±2z/ + 2n) 



(73) 



2(2-e) 5^ X --±u;-±u; C^j + ((2-e)7rtan(7rz/) + rf^jC^ 



where in the last step we chose the integration constant to be 

;i-60'(2- 



cte(z/) 



647r2 



(2 — e)7r tan(7rz/) + d^ 



(74) 



(2) 

where the tangent contribution is fixed by requiring that pg be finite for all z/, and a 
u dependent function, which is finite for all u. Notice that the divergences at half integer 
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V coming from the first hypergeometric function are cancelled by the contribution from the 
tangent in (1711) . The integration constant (I74p scales as radiation, oc l/a^, and it can 
be uniquely fixed by calculating the full stress energy tensor instead of just the trace. But 
in order to do that, we need the propagator off coincidence, which we do not have at this 
moment. 

With this we are now ready to calculate the quantum energy density and pressure, which 
are obtained by integrating Eq. fl5B]) with the help of Eqs. (11371) in Appendix B. The result 
for the energy density is 



,m = (l-60^(2-e)/7^ 



1, / Attu^ \ Cy 6-7e 
m 



2 V(l-e)^^^(l + C)V 2 12(1- 
± 



- E X ^'■i (i. i ± i ± q + 2 

-2(2-£)C-'l.i(ii^)+4(l-£,^^^ 

(-■i±'lv /II \ 

-2(2-^)E^x=^'.(i.-5±''^5±''^C') 

+ ((2-e)7rtan(7rz/) + d.jC' - 4(l-e) E IZ^j 



(75) 

with 1/ > 0. The total energy density is obtained by adding this to the ultraviolet contribu- 
tion (l69l) . The quantum pressure is then simply a sum of 

= ^(Tf +pf) (76) 

and the ultraviolet contribution y^q^ in Eq. (!70l) . The equation of state parameter Wq is , 



Wq = ^. (77) 

In other to find out whether the quantum contribution can become important, we need to 
investigate whether the criterion fl6T|) for the growth of quantum contribution with respect 
to the classical contibution is ever met (see also Eq. fl62|) above). Rather then studying Wq{t) 
in its full generality, we shall consider only the leading (late time) behaviour of Wq. For 
this we need the leading order behavior of Tq, pq and Pq. Remember that in this section we 
consider an accelerating universe such that at late times ( ^ 0. 
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7.1.1 The case when e < 1, < u < 3/2 



From Eqs. ([M]), ([ZI]), ([69]) ([70]), ([75]) and ([7^ follows that when < z/ < 3/2 the leading 
order late time contribution to the one loop stress energy trace, energy density and pressure 
are, 



Pi 



(l-60^(l-6)(2-6)g^ 

;i-60'(2-e)ij'^ 

:i-60^(2-e)g^ 
647r2 



6e 



6e In ( ^ 



In 



Ho 
H 



4fl- 



6e 



(78) 
(79) 
(80) 



where 



HHo) 



In 



47r/i^ 



3-5e 



TC 



2-e / 1 



2 6(l-e) 
1-e 



4(l-e) 
11527r2 



tan(7rz/) 



+ 



3(l-e)V4e (2-e)(l-60/ (1-60^ 
These relations then imply for the equation of state parameter Wg = Pql Pq- 



^ 3 / 3ln {Ho/H) ^ ^ 3^ 



(82) 



The last implication follows from the observation that at late times a —>■ oo, H oc a^*^ ^ 0, 
such that formally ln{Ho/H) ^ oo as long as e > 0. When e = one recovers the well 
known result, Wq = —1. The results fl78] - [82|) also apply when u is imaginary, or more 
generally whenever ^[u] < 3/2. 

7.1.2 The case when e < 1, u > 3/2 

When u > 3/2 in general the lowest u dependent power in (ITTl) and (1751) dominates the stress 
energy tensor at late times. From the more convenient form (17214731) expressed as series we 
can read off the dominant contribution to Tq, pq and Pq-. 
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-J- n 



167r2 
(l-60'(2-e)if^ 
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647r2 
647r2 



6e 



3-2z/ 

6e 
3-2z/ 



3-2u 

(3-5e 



+ (3-5e) - (l-e)(3-2z/) 

4C^ 



(l-e)(3-2z/ 
(3-5e) - (l-e)(3-2z/ 



c 



3-2i/ 



l + 2u 

4[e-2(l-e)z/; 
3[l + 2z/l 



(83) 
(84) 
■C=^-2-(85) 
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From this we conclude, 



Wg^-^-\{l-e)u+'- (3/2 <u) 

(86) 

^ + 3 ((3 - e) - v/(3 - e)2 - 24(2 - e)^) . 

from which we see that Wg < Wf, if ^ is negative and Wg > Wbis ^ is positive, for all relevant e 
(remember that this analysis holds for e < 1). Thus we find that the quantum contribution 
to the energy density grows with respect to the backgound energy density for ,^ < 0. 

This result is correct however, provided the leading order contributions (I83l485|) do not 
vanish, which is the case when the expression in the square brackets does not vanish, i.e. 
when 

6e 

__ + (3_5e) - (i_e)(3-2z/) ^ . (87) 
We shall now study in some detail the special case when flHTl) vanishes. 



7.1.3 The special case when e < 1, ^ = 

It is instructive to observe that ( 1871) vanishes for the following two values of u: 

3-e 



or equivalently when 



\ (3- 



-e 



i2 



The first value in (IHHl) is irrelevant since we consider v > 3/2. The second value however 
might be interesting. It corresponds to the minimal coupling of the scalar field to curvature, 
.^ = 0. If3/2<z/<5/2 (or equivalently < e < 1/2), the leading order contribution in this 
case is logarithmic and the above analysis presented in subsection 17.1.11 applies and we have 
(c/. Eq. 



Wg^-l + -e (^ = 0, 3/2 < z/ < 5/2, < e < 1/2) . (90) 
3 

When, on the other hand, z/>5/2(l/2<e<l), the dominant contribution comes from the 
subdominant term in the u dependent series f l72H751) . Since we study the case where ^ = 



4 
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we have from fll8l) that v = (3 — e)/[2(l — e)] and we can write the stress energy contributions 

as, 



r C-^o (l-60^(l-6)^(2-6)(l+6)/7V s-2. .Q^^ 
" 167r2(l-2e) ^ ' ^ ^ 

^ - ^-''^^St^"'^''^ ^'-'' (^/^ < - < - ^) ^^^) 

^ (l-60^(l-6)^(2-6)(l + 6)g^ 

967r2(l-2e) ^ ' ^ ^ 

Thus the equation of state parameter is in this case, 

w;,^-^ (e = 0,z/>5/2,l/2<e<l). (94) 

In figures [3l H] and [5] we illustrate the results of this section (these figures also already 
show the results for the decelerating case, e > 1 which we shall discuss in section [7]2]). From 
the previous discussion we found that in an accelerating universe quantum effects can become 
important if < while they will always be subdominant for ^ > 0. In figure [3] we show the 
borderline behavior: ^ = 0. The equation of state parameter of quantum fluid Wg = — l+4e/3 
(when < e < 1/2, conform (pOl) ) and Wg = —1/3 (when 1/2 < e < 1, conform (1941) ) is 
always larger than the corresponding background parameter Wf, = — 1 + 2e/3, implying that 
the criterion (l6Ti) is never met and the quantum contribution can never become important. 
If, on the other hand, ^ 7^ the situation changes drastically. In figure H] we show Wg for 
positive values of ^. The curves are given by = — 1 + |e for u < 3/2 (conform fl82|) ) and 
Wg = — 1(1 — e)z/ + |forz/ > 3/2 (conform fl86l) ). We indeed see that we always have Wg > Wb, 
such that the energy density due to the quantum effects dilutes faster then the background 
energy density. In figure [5] we show Wg for negative values of C,- Since for negative ^ we 
always have that > 3/2 we have (conform fl5^ ) that w^ = — 1(1 — e)z/ + |. We indeed 
see that Wg < in this case and thus the energy density in the quantum contribution will 
dominate at late enough times. The time it takes for the quantum contribution to become 
important can be estimated from Eq. (l62l) to be 

1 ( Mp\ -6?(2-£) 



—^p;^ ^ (e<0,0<e<l). (95) 

For < —6.^ ^ 1 this time can be very long, much longer than the age of the Universe. In 
fact —6.^ < 1 can be tuned to make tf. of the order of the age of the Universe, in which case 
one could relate the contribution of quantum fluctuations to the dark energy of the Universe, 
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and perhaps even address the question 'Why now!' of the dark energy dominance. The above 
analysis works however only for matching radiation era onto accelerating Universes, in which 
the Universe spends most of the time, and that does not correspond to the observed Universe. 
Therefore, we also need to consider radiation matching onto the decelerated universes, which 
is what we hurriedly do next. 



7.2 Matching onto deceleration for a general u 

The exact result for Tq'^^ in Eq. (!7T!) is also correct for the decelerating universe. Since in a 

(2) 

decelerating universe at late times C — > oo, it is more suitable to express Tq in terms of the 
hypergeometric function of The result is. 
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(96) 



2(1- 



where we made use of the identity, 

^2a 



2a 



X2F1 l,a;a+l;C 



/•2{q-1) I 



)+^cot(7ra) (C'>1). (97) 



In the derivation of this identity, we took the central value along the cut of the hypergeometric 
function ("^ > 1 (this prescription removes a possible imaginary contribution). Notice that 
for a = (3/2) ±z/ the last bit in (^71) does not contribute, since J2± cot[7r(3/2) ivrz/] = 0. The 
result fl96|) can be also obtained from Eq. fl5^ by acting with □ on Eq. fj471) . representing a 
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check of ( !96|) . The identity fl97j) can be also used to derive pf^ from Eq. ([75 
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The one loop pressure pq is, as before, obtained by inserting ( |96l) and (|98l) into relation (|76l) . 
We are now ready to consider the late time limit, a ^ 00 which implies for the decelerating 
case that Q ^ 00, oi the stress energy tensor and the corresponding equation of state 
parameter in a decelerating universe. Before we begin analysing particular cases, let us 
observe the general structure of Eqs. and flUS]) . The curly brackets in contain terms 
that grow oc and terms with z/-dependent powers of The z/-dependent powers can be 
summarised as. 
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Similarly, pf^ in ( |98l) contains terms that grow as oc C^, and the terms that contain v- 
dependent powers of C, which can be summarized as. 



[l-Qif{2-e)H' 
64^^2 



EE 



± n=0 



6e 



-l±2z/+2n 



(3-5e) - (l-e)(-l±2z/+2n) 



4C 



l=F2;/-2n 



3 ± 2z/ + 2n 
(100) 

This form could be also obtained by integrating fl99|) . representing a nontrivial check of 
our result fl98|) . It is now clear that the leading order late time behavior of pq and Tq 
is dictated either by the term oc C,^ in pq or by a term oc (^1+22^^2" ^j^]-^ _ 0, 1, . . . in 
Eqs. f l99] - ll00p . In both cases the ultraviolet contributions fl68]470|) are subdominant and can 
be neglected. Let us consider first in somewhat more detail the contribution oc C^. This 
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contribution is the dominant one, for all z/ < 3/2. It can be easily seen to lead to a = 1/3 
(since this contribution does is not influenced by the trace of the stress energy tensor, its 
contribution must be traceless), independent of e, and therefore this contribution will always 
lead to a strong backreaction, for larger enough values of e (such that Wb is large). We 
do not believe however that this effect is due to infrared particle production. The reason 
is the following. A large backreaction, on physical grounds, is only expected when due to 
particle creation the infrared becomes highly correlated. This high amount of correlation 
leads to the infrared divergence and thus is present only for u > 3/2. Now not only does our 
procedure not fix this contribution uniquely, in Appendix C we show that this contribution 
cannot be determined uniquely in the present model. The reason is that it turns out that 
the contribution to the stress-energy tensor oc H'^('^, is ultraviolet divergent. The physical 
origin for this divergence is probably the sudden matching between the two space-times 
we consider. Such a divergence then needs to be renormalized, and this renormalization 
inevitably introduces an arbitrariness in the coefficient in front of the H'^C^ term. Therefore 
we are free to choose the constant d^, and for the purpose of present paper we choose the 
constant di, such that the term proportional to cancels. In this case we have that when 
u < 1/2 a term proportional to ("^ dominates, or when u > 1/2 the dominant term is 
proportional to (^i+2!^-2n._ 



7.2.1 The case when e > 1, < z/ < 1/2 



In this case the dominant contribution to the one loop stress energy tensor comes from the 
term oc ("^ in (p6!) and 



Pg 
Pg 



(l-60^(l-6)(2-e)g^ 

' 167r2 
C-^ (l-6e)^(2-e)g^ 

' 647r2 
C-«o {l-60\2-e)HU 2 

' 647r2 



l-2e) 



from which we conclude: 



>oo 



1 2 



Wb + ->Wb, 



(101) 
(102) 
(103) 

(104) 



where Wb denotes the background equation of state parameter ( 158|) . Notice that the condition, 
z/ < 1/2 implies that ^<l/6, ife>2 and ^ > 1/6, if 1 < e < 2. In those regimes we 
therefore do not expect any significant backreaction. 
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7.2.2 The case when e > 1, z/ > 1/2 



In this case the dominant contribution to Tg, pg and Pg comes from a u dependent power in 
Eqs. ([96]) and ([98]) (c/. also Eqs. MTOOD): 



T 

■J- n 



Pg 



Pq 



(J— >oo 



(l-60^(l-6)(2-e)ff^ 



167r2 

C-K^ (1-60^(2-6)^^ 
' 647r2 



6e 



6e 



(5e-3) + (e-l)(2z/+l) 



f — >oo 



647r2 



l + 2z/ 

6e 
l + 2z/ 



l + 2z/ 

- (5e-3) + (e-l)(2z/+l 



C 



1+2^ 



(5e-3) + (e-l)(2z/+l) 



3 
4 

3 L 



- 2u 
fe-11 



(105) 
(106) 



1 



3-2Z/J 



(107) 

where we assumed that the term in square brackets does not vanish and we took u ^ 3/2 [f]. 
From Eqs. (I105H107]1 we conclude: 



2 2 

---(e-l)z/ 

3 3^ ' 



-(e-l)z/, 
3^ ' 



(108) 



This equation implies that the quantum contribution will eventually dominate provided 

e + 1 



V > 



2(e-l) 



(109) 



or equivalently when 



e < 

e > 



e-l 



3(2-e) 
e-l 
3(e-2) 



when 1 < e < 2 
when e > 2 . 



(110) 



Together with the criterion ^ < when < e < 1, these relations define the regions 
plotted in figure [6l for which one expects that the quantum one loop contribution to the 
stress energy will eventually dominates over the classical stress energy tensor that drives 
the Universe's expansion, and in which case the quantum backreaction can infiuence - and 
in fact change - the evolution of the Universe. These conclusions are correct, provided the 
contributions fll05H107|) do not vanish, which will be the case when the sum of the terms in 
square brackets does not vanish, which is the case we consider next. 



•^Whcn v = 3/2, pq and Pq acquire a logarihmic contribution of C: which yields the same Wq as is in 
equation (I108|) . 
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7.2.3 The special case when e > 1 and ^ = -(e-l)/[3(2-e)] 

The case when the sum of the terms in square brackets (I105H107|1 vanishes requires a special 
attention, which is the case when, 

\ '2(e-l)j ^ 24(2-e) ' 3(2-e) 

We shall only consider the second value, since we do not expect anything significant for 
u < 3/2. For this case the results of the subsection 17.2.11 apply for all < < 3/2, and we 
have. 




12 2 

Wg -- + -e = Wb+->Wb (0 < z/ < 3/2 ^ e > 2) , (112) 

When, on the other hand, 1 < e < 2, such that u > 3/2 and ^ < 0, the subleading terms in 
the sums fl99l - ll00l) dominate, 

y (l-60^(l-6)^(2-6)(5-36)g^ ^,^(,_,^ 

^ C-oo (l-60^(l-6)(2-6)(5-3e)g^2(e-l) ^,/(,_i^ 

647r2 2e-3 ^ ^ ^ 

(l-60^(l-6)(2-e)(5-36)g^ 2(e-l) 
' 647r2 3(2e-3) 

This then implies for the equation of state parameter, 

„,._|H-i.^.. + |(.-l) (.>3/2«l<.<2). (116) 
Together with Eq. flll6p this implies that Wq > Wb, Ve > 1. 

Let us now go back to the figures [3], H] and [5l already discussed at the end of section 17.11 
Now we can also understand the part of the plots for e > 1. For the decelerating case we 



2(2e-3) + l (2/(^-1) _ (115) 



found that quantum contributions can dominate if ^ < —3^1^ for 1 < e < 2 and ^ > 3(2-^) 
for e > 2. In figure [3] we show the borderline case, thus for e > 1 we choose ^ = — • 
This implies that for 1 < e < 2 we have that = — | + |e, conform (11161) and for e > 2 
we have ly, = — | + |e, conform (11121) . In figure H] we show Wq for positive values of ^. A 
positive value of ^, less then 1/6 means that u > 1/2 for e < 2 and u > 1/2 for e > 2 thus 
in this case we show Wq = e — | — |(e — l)z/, conform (11041) . for the region 1 < e < 2 and 
Wq = + |e, confrom (llOSp for e > 2. On the other hand if > 1/6, we have that for 
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1 < e < 2 that z/ < 1/2 and for e > 2, i/ > 1/2. Thus in this case it is precisely the other 
way around. Moreover , since now ^ > ^~^2) region where e > 2, we find conform 

f ill 01) that in this regime Wq can become smaller then Wb. Finally in figure we show Wg 
for negative values of ^. If ^ is negative, we always have that u > 1/2 for 1 < e < 2 and 
z/ < 1/2 for e > 2. And thus we show Wg = e — | — |(e — l)z/, conform (11041) . for the region 
1 < e < 2 and = — | + |e, confrom (11081) for e > 2. We find, conform (IllOp that in the 
region 1 < e < 2 we can have that Wg < Wb- 

This completes our discussion for now. However there are special points, corresponding to 
u = 3/2 and u = 5/2, where the above discussion fails because of the logarithms appearing 
in pq and Pg. However it turns that the scaling of the quantum corrections, Wg does not 
deviate from the behavior described above. Details on this we present in Appendix C. 



8 Comparison with the cut-off regulated T, 



In earlier work [18] we have calculated the one loop stress energy due to scalar field fiuctu- 
ations, where the infrared was regulated, by working on a compact spatial manifold, or in 
other words, by assuming that the Universe is a comoving box, with periodic boundary con- 
ditions (which to a good approximation can be described by a comoving infrared momentum 
cut-off). 

Generalising the result from [18] to a nonminimally coupled scalar, we obtain for the 
renormalised stress-energy tensor 



(0|T^.|0)=- 



6(2-6)(l-60^ 
IGvr^ 



l-2u' 1 + 2u' 

+ 2 + 2 

l-2v l + 2v 



lE + In 



+ ^(^-z/) + ^(^ + z/) 



167r2 



1 - 40(2 - 6)6^ - ^(1 - 60(1 - 4e + ^)e]aHyu 



+ ( - i( - 7 + 8e(l - 40 + 30e) - ^"^(3 - 22e + 22t^ - At^))g^, 



(117) 



where v' = dv / dD\£,=i = and the terms in the last line are the infrared corrections arising 
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{Wq, Wb) 




Figure 3: Wq (blue) and Wb (red, dashed) versus e. In this plot we choose ^ to be on the 
boundary value which seperates regimes where Wq can become less then Wb and regimes where 
this can not happen. This implies ^ = for e < 1 and ^ = — ^'^^ e > 2. Wb is given 

by Wb = —I + |e and Wg is given by = — 1 + |e (0 < e < 1/2), Wq = — | (1/2 < e < 1), 
w,; = -| + |e (1 < e < 2) and Wq = + |e (e > 2). 
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Figure 4: Wg and Wh (red, dashed), versus e. This plot shows positive values of ^: ^ = 0.1 
(blue, solid) and C, = I (green, dotted). The plotted values for Wg are for the ^ = 0.1 curve: 
Wg = —1 + ^e, for < e < ec, where ec is determined by the requirement u = 3/2, implying 
ec = 0.44 Wg = -\v +1 (ec < e < 1), = e - |(1 + z/), (1 < e < 2) and Wq = -\ + |e 
(e > 2). For the case where ^ = 1 we have = — 1 + |e (0 < e < 1), Wg = — | + -^e 
(1 < e < 2) and = e - |(1 + v), (e > 2) 
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{Wg, Wh) 



2 



-1 




-2l" 

Figure 5: Wg and Wb (red, dashed), versus e. This plot shows negative values of ^: ^ = — 0.1 
(blue, solid) and ^ = —0.5 (green, dotted). The plotted values for Wg are for both curves 
Wg = -|z/ + I (0 < e < 1), = e - |(1 + z/), (1 < e < 2) and = -| + |e (e > 2). 
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Figure 6: Boundary values of as a function of e. The solid blue regions are the regions for 
which the quantum stress energy scales slower that the classical stress energy, such that the 
quantum contribution will eventually dominate over the classical contribution. The regions 
are bounded by ^ = (0 < e < 1) and ^ = — ^ ■'-)■ '^^^ dotted asymptotes are given 

by e = 2 and ^ = 1/3. 
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from the infrared cutoff, which can be written as [IS] 



(0|T^i.|0)Ar = (pn +PN)a KK+PnGixv 



1 (2(1 -e)(Ar-z/) +3-e) T{u ~ N)T{2v - N) 



4vr5/2 

X hHi 
'1 

X 



3 + 2N -2v 



60 (iV- 



V 



A^)r(A^ + 1^ 



4 



3Vl-2iV + 2i/ 



(1 



1-2N + 2U 



(118) 



and 



(0|T^.|0) 



N 



1 (2(1 - e)(Ar + I/) + 3 - e) r(-z/ - iV)r(-2z/ - AT) 



47]-5/2 

X i/^(l 
^1 



3 + 2Ar + 2z/ r(i-z/ 

3 VI 



X 



A^)r(A^ 



3\1-2N -2u 



1 



2N -2u 



(119) 



where Zq = kolrjl, is the (comoving) infrared cut-off scale (fco = 27r/L, where L is the 
comoving size of the Universe) and t] is conformal time, defined by ar] = — 1/[(1 — e)H]. 
We assume that the size of the Universe is at an initial time tjq super-Hubble, such that 
^ol^ol = ko/[{l — e)aQHo] <C 1. Notice that in an accelerating universe (0 < e < 1) at late 
times zq decreases, such that — ^ as ?7 — > 0— (or equivalently a oo). On the other 
hand, in a decelerating universe t] increases, resulting in an increasing zq. When zq ^ 1, 
the size of the Universe becomes comparable to the Hubble radius. Even later the comoving 
box shrinks to sub-Hubble sizes, and Zq ^ oo when rj^a ^ oo. This case requires a special 
attention, and it is treated below. 

We shall now construct the leading order contribution from the corrections flll8Hll9l) to 
the stress energy tensor in both accelerating and decelerating universes. If Zq approaches 
zero, we see that the leading order contribution comes from the = term of (11181) . This 
term is growing if i/ > 3/2, as is expected, since this is the requirement for an infrared 
divergence ( l25i) . The case v = 3/2 leads to a logarithmic growth which, for brevity, we do 
not study here. The results turn out to be analogous to the case we consider here, described 
in Eqs. (I140H142"]) for the accelerating case. If v > 3/2, the leading contribution in an 
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accelerating universe is 



(0|T,..|0) 



4:H\1 - ey 



37r2(l + 2z/)(2i/-3) V2/ 



(-) (i-60r(z/)r(z/ + i)^ 



3 + e 
2(1 - e 



2(1 
+ 0(^, 



(120) 



J 



This expression vanishes when ^ = 0, since then z/ = (3 — e)/[2(l — e)]. In that specific case 
we need the next to leading order contribution, which is proportional to Zq~'^'^. The form of 
the leading order contribution (I120p very similar to corresponding matching case (!83] - [86|) . 
with zq ^ ( (the precise numerical coefficients multiplying the leading order terms are, of 
course, different). Thus in an accelerating universe, the two regularization procedures are 
qualitatively the same. The reasons for this is that any infrared regularization implies that 
we effectively suppress modes with wavelengths larger then some scale, be it given by zq or (. 
In an accelerating universe, this scale grows faster then the Hubble radius and therefore the 
precise details of this regularization become less and less visible as time goes on. Therefore 
we indeed find that at late enough times the two regularization schemes give qualitatively 
the same result. 

However, as mentioned above, in a decelerating space-time we have a different situation: 
the sums over N run to infinity and, since zo grows, this leads in principle to fast growing 
terms. Moreover we see that this happens when the Universe's size becomes sub-Hubble, 
independent of z/. Thus also infrared perfectly finite space-times become dominated by the 
cut-off. The most reasonable approach is is this case to sum the sums over N when 2:0 -C 1 
and then analytically extend to late times when 2:0 S> 1 and the Universe's size is sub-Hubble. 
We shall now illustrate how this procedure works by calculating the trace of the stress energy 
tensor. After taking the trace of (11181) we can perform the sum over to obtain 



N 



N=0 



.^(5-2.)(3-2^(l-«««^(l-^''(f) 



3-2!/ 



i']'2F^{ - - u, - - u-A ~ 2^, - - -z/; -z, 



' '35.7 



») 



+ ^{3 - 2z/)2F3 (- - z/, - - z/; 2 - 2z/, - - z/, 1 - z/; -z, 



'1211 



The leading order can be studied by considering the asymptotic expansion of the 2F3 hyper- 
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geometric functions. We have in general 

Tih)T{b2)Tih) 



2F3 01,02; 61, 62,^3; -z 



T(ai)T(a2-ai) 



+ 



+ 



T{a2)T{ai - 02) 



T{bi - 02)1(62 - a2)r(63 - 0-2) 



r(ai)r(a2) r(6i-ai)r(62-ai)r(63-ai) 



(l + 0{z-^) 



cos (nx + I^Tz) + ((3ai + 3a2 + 61 + 62 + 63 - 2)(8x - 2) 

+ 16(6162 + &1&3 + &2&3 - 01^2) - 3) sin (txx + 2v^) (1 + 0{z-^) 



with 



X 



^ {ax + 02 - 61 - 62 - 63 + 



(122) 
(123) 



Using this we find that the leading order terms are 



(l-60(l-e)3ij-^ 
Stt^ sin^fvri^) 



Ar=0 



^0 , 2/5e-l 



1 — e)z/^) sin (22:0 + 7rz/) 



— Zg(l — e) cos(2zo + 7rz^) 



(124) 



The second series Y^=q^\^^ tX^)^ -^l- (Hl^p can be obtained from fll24p by interchanging 
V with —V and can be easily added, resulting in. 



V / Stt"' sm tti^ 



Af=0 



2 _|_ 2 
+ 



5e-l 



+ 2(1 — e)z/^ ) cos(7rz/) sin(22;o) — 2zq(1— e) cos(7rz/) cos(22;o) 



Selecting the leading order term in (11251) we can recast it as. 



T ~ 

J-q — 



;i-60(e-l)gofcg cos(7rz/) 
471^ sin^fvrz/) 



cos 



2A;n 



:e-l)ifa 



-e-3 



(125) 
(126) 



We shall not attempt to evaluate Pg and in this case, since the asymptotic expansion would 
result in an expression that depends on the lowest value of zq. Instead we shall compare 
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the corresponding traces. When Eq. fll26p is compared with the background contribution, 
T^ = p^^- Sp^ oc l/a3(i+'"'') = one finds that 



which grows when e > 3. Thus we find that the hmiting case is kination, for which e = 3 
and pb oc 1/a^, and the quantum contribution to the trace Tg given in (11261) scales the same 
as the background contribution to the trace. For all e > 3 we thus find that the quantum 
contribution will eventually dominate over the background energy density. There is one 
exception: as long as u is not half integer, the scaling (I127P is correct. However when u 
is half integer, cos(7ri/) = and the leading contribution (I126p vanishes. The dominant 
contribution is then the term oc in Eq. (I125p . We have seen in section [7.2. II that such a 
contribution will never lead to a strong backreaction. This is in contrast with what we found 
in the present work. Provided that we correctly disregarded the term, which arose as an 
integration constant, we have found in section [7^ that, if the infrared is regulated using the 
mode matching, for all u > 1/2, the scaling of the quantum energy is indeed governed by u 
and not e. 

9 Summary and discussion 

In order to facilitate the reading of this rather technical paper, we shall now recap our main 
results. The main motivation for this paper is to study the role of the quantum infrared 
fiuctuations generated by the expansion of the Universe for massless scalar fields, with a 
possible coupling to the Ricci scalar. 

It is well known that the infrared (IR) sector of such a scalar field posses problems on 
a cosmological background: in space times with a negative pressure, and constant accel- 
eration/deceleration parameter e = —H/H"^, the Bunch-Davies (BD) vacuum of massless 
minimally coupled scalars is infrared singular [4j (see section [2]). 

Of course, we know that our Universe must be infrared finite. That means that the 
infrared sector of the theory must be regulated. It is currently unknown what is the precise 
nature of the regularisation scheme realized in our Universe. Hence it is worth investigating 
different regularisation schemes and compare the results. In the end however, the presence of 
this infrared singularity is precisely what makes the present study interesting. The singularity 
indicates that there is a growth in long range correlations due to particle production. Even 
after regulating the infrared divergence, this growth of the correlations is still there, since 
it is simply a physical effect. These growing correlations can then in principle - after a 




(127) 



40 



sufficient amount of time - contribute significantly to the energy density in the Universe. If 
this is so, they might change the evolution of the Universe significantly. 

In this work we focus our attention on the regularisation of the infrared which is executed 
by an epoch of the very early universe which has very little particle production and whose 
BD vacuum is therefore infrared finite. For definiteness we choose this early epoch to be 
radiation era. We then match it onto a constant e space time, and calculate the corresponding 
coincident scalar propagator. The exact solution for the mode functions can be expressed 
in terms of Hankel functions (fT9l) with a rather complicated index v (fT8|) . The BD vacuum 
corresponds to the choice (122|) . implying that one considers only positive frequency modes. 
As can be seen from Eq. ( !25|) the propagator for the BD vacuum is IR singular whenever 
ly > {D — l)/2, where D denotes the number of spacetime dimensions. At the matching 
we require continuity of both the mode functions and their ffist derivatives. This implies 
that after the matching, the mode functions become a mixture of positive and negative 
frequency modes, in such way that the (coincident) propagator fl32l) is IR finite for all e. 
The propagator still suffers from the standard (logarithmic) ultraviolet divergence, which in 
dimensional regularisation appears as a 1/(D — 4) divergence (l36l) . This divergence induces 
an analogous divergence in the stress energy tensor flM|) . and can be removed using standard 
techniques by an counterterm flM]) . 

To construct the propagator, we have ffist calculated the coincident propagator for half 
integer u (140] - HT|) and then analytically extended the result to all (complex) z/ fHSl) . Since 
the coincident propagator for half integers (142!) is valid for all times later than the matching, 
we suspect that the analytic extension (11^ is unique. From the coincident propagator, one 
can derive the trace of the one loop stress energy tensor Tg, based on equation (IMIl . The 
symmetries of the background space time dictate the perfect fluid form of the quantum 
contribution to the stress energy tensor and using the covariant conservation of the stress 
energy tensor fl56|) . one can derive the individual components (quantum energy density 
Pg and pressure Pg) of the one loop (T^u)q- The general results of this procedure are for 
accelerating universe presented in Eqs. flTTl) and in Eqs. flTSHT^ and for decelerating universes 
in Eqs. ( !96] - [98|) . Since we obtain pg and Pg by an integration procedure, our results are 
unique up to an integration constant. Adding this integration constant corresponds to adding 
a component to pg and Pg, which scales as a radiation fluid, oc and thus does not 

contribute to the trace of the stress energy tensor Tg. To fix this component uniquely, we 
would need to know the propagator away from the coincidence limit, which we have not 
derived in this paper. Since this undetermined radiation component yields a subdominant 
contribution at late times for accelerating universes and thus it is of no relevance for our 
discussion of the late time quantum stress energy tensor in section 17.11 It may be however 
relevant in decelerating universes. For the purpose of this work we have ignored this radiation 
contribution in our analysis of the late time quantum stress energy tensor in decelerating 



41 



universes presented in section 17.21 

In order to study the significance of tlie quantum one loop stress energy tensor, in sub- 
sections [7?T] and [7]2] we study in detail tlie equation of state parameter Wg = Pq/pq both in 
accelerating and in decelerating universes. If the following criterion is satisfied fl6Tl) . 

Wq < Wb, 

(where Wb = Pb/pb is the equation of state parameter of the background fluid, driving the 
expansion of the universe), then the quantum contribution to the energy density dominates 
over the background contribution at late times. In this work we do not attempt to analyze 
what exactly happens in such a case (to start with, the assumptions underlying our calcula- 
tions would become incorrect), but simply want to answer the question if and under what 
conditions this criterion is ever satisfied. 

Our results are presented in detail in sections 17.11 and I7.2[ Since several cases require a 
separate discussion, it is difficult to get a quick grasp of the results. In order to facili- 
tate a quicker understanding of our results we present in figures |3HS1 the quantum equation 
of state parameter Wg vs the background equation of state parameter Wb = —I + (2/3)e. 
The criterion fl6T|) then tells us when the quantum contribution becomes dominant over the 
background contribution at late times. We shall now describe our results in some detail. 

Firstly, in figure [6] the shaded regions represent the regions in parameter space {^, e} 
where the criterion Wq < Wb is satisfied. For accelerated universes (0 < e < 1) this means 
simply that when ^ < 0, Wq < Wb- For decelerating universes (e > 1) we need to distinguish 
two cases. When 1 < e < 2, then ^ < — (e — 1)/[3(2 — e)] flllOp assures Wg < Wb- When on the 
other hand e > 2, ^ > (e — l)/[3(e — 2)] (II 101) is required in order that Wg < Wb- Keeping in 
mind that the Ricci scalar curvature R = 6{2 — e)H'^ is positive when < e < 2 and negative 
when e > 2, we see that Wq < Wb can be met only when the effective scalar 'mass' parameter 
m^jj = is (sufficiently) negative. A careful analysis presented in subsections l7.1.3l and l7.2.2l 
shows that at the boundary of the shaded region in all cases Wg > Wb is satisfied, such that 
the quantum contribution can never become important. This can also be seen from figure [31 
In figure m we show Wq vs Wb for positive ^ > 0. When the coupling ^ is smaller than 1/3, 
we have that Wq > Wb such that the quantum stress energy can never be important. When 
however ^ > 1/3, in spacetimes with e > 2, Wq can become smaller than Wb, indicating the 
late time dominance of the one loop contribution. 

Finally, if ^ is negative, which is shown in figure [51 In this case Wq < Wb is always met for 
accelerating cases, but also for decelerating cases whenever ^ < ^cr = — (e — l)/[3(2 — e)]. 
Notice that the critical e (for which ^ = ^cr) is always smaller than 2. 

Apart from the regularisation scheme presented in this work that involves matching 
from a nonsingular spacetime, other infrared regularisation schemes are possible. We do 
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not know which regularisation scheme was chosen by our Universe. We do know however 
that, since infrared regularisation is physical, there are in principle physical observables that 
can distinguish between different regularisation schemes, and thus we should be able (at 
least in principle) decide which one was picked by our Universe. But in order to find out the 
answer to that question, we need to investigate different plausible IR regularlisation schemes. 
Other regularisation schemes include: placing the Universe in a large (comoving) box (this 
corresponds to an infrared momentum cutoff) which has been explored in Refs. [TT t [T51I52] . 
a tiny scalar mass (it is not clear whether that is possible to implement for the graviton), 
a positive spatial curvature, subtracting an adiabatic [20j or a comoving vacuum [221 l2T] . 
etc. While most of these schemes are physically well motivated, the implementation is often 
hindered by our lack of knowledge of the relevant propagators. In particular, we know the 
propagators in positively curved universes and for massive fields only in very special cases 
(de Sitter space, radiation era), which is not enough to conduct a sufficiently general analysis 
of the quantum backreaction regulated this way. 

Owing to the fact that the analysis of the Universe in a finite comoving box can be well 
approximated by an infrared momentum cutoff, in Ref. [18] we were able to perform the one 
loop analysis of the massless scalar backreaction in expanding universes with a constant e 
parameter. In section [H] we present in some detail a comparison between the two regular- 
isation schemes. We find that in accelerated universes the leading order contributions to 
the corresponding late time stress energy tensors are of a similar form in both regularisation 
schemes in the sense that the leading order behavior with the scale factor is identical (c/. 
Eqs. ( I55H5BI) and (11201) ). albeit the coefficients of the leading order terms differ. In fact, the 
coefficients are comparable when the cut-off scale is chosen to be equal to the horizon scale 
H at the matching, i.e. when the initial comoving size of the Universe is of the order the 
Hubble radius. This is understandable, given the fact that the comoving box in accelerating 
universes grows with respect to the Hubble radius, i.e. as times goes on the Universe's size 
becomes more and more super-Hubble. The precise details about what precisely happens on 
super-Hubble scales become less and less 'visible' as time goes on. This is precisely what we 
find: at late times the two approaches give qualitatively the same answer. 

On the other hand, in decelerating space-times the matching and cut-off regularisation 
schemes yield very different results. The leading order contributions to the trace of the 



stress energy tensor are in this case (c/. Eqs. (llOSp for u > 1/2, Eq. (1 10 II) for 3?[z/] < 1/2 

and ([nsHns])) 





(matching, 3fJ[z/] < 1/2), 



(cut — off) 



(129) 



(128) 
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where ( = aH/{aH) and zq = k^/Ke — l)aH] = {k^/Ke — l)aH]}(, such that Zq/C is a 
constant given by the ratio of the physical cutoff scale ko/d and the (inverse) particle horizon 
~ (e — 1)H at the matching time. For the cut-off regulated case we find that the growth is 
independent of i'. The quantum contribution in this case becomes more and more dominant 
if e becomes larger. In the mode-matching case we see that the growth is dependent on u 
and moreover the effect is more profound the greater u is. This is exactly what one would 
expect based on fl^ . What happens physically is that, in a decelerating space-time, physical 
scales grow slower then the Hubble radius, and eventually the comoving size of the Universe 
becomes sub-Hubble. This would eventually dominate all effects and thus what one is looking 
at then has nothing to do anymore with infrared particle production! The results (11291) 
and (I125H126]) are obtained by analytically extending to the limit Zq ^ oo, which is precisely 
the limit in which the Universe's size becomes sub-Hubble. Thence, not surprisingly, in this 
limit Tq oc 1/Vc becomes a function of the comoving volume Vc of space-time, representing 
an observable indicating the size of the Universe. In the infrared reularization presented in 
this work, all potentially relevant effects are however due to infrared particle production, 
and since this is the physical effect we wish to study, we feel that this approach therefore is 
advantageous in decelerating spacetimes. 
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Appendix A: Acting with the d'Alembertian 

In order to calculate the contribution to the trace of the stress energy tensor (I54p . the 
following identities are useful, 

mH'fiO] = (l-^)^'[6^+(3-5e)C^-(l-e)C^C^]/(C) (130) 
n[HHn{H^)] = 2eH^\3{l-e)\n{H^) + {3~5e)], (131) 
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where we made use of □ = —{dt + SH)dt and dt = —{l — e)H(di^. From these we easily get 
the following useful identities, 



□ [if^ln(l±C)] 

n[HHn{C)] 



:i-e)H' 



6eln(l±C) + (3-5e) - 



2(2-3e) 



1-e 



i±C (i±C)^ 



6eln(C) + (3-5e 
6e + (3-5e)cj- (l-e)cu2 



(132) 



From these identities we can also obtain how the d'Alembertian acts on the logarithms 
in (IM 



1\2 



'l-e)H''{ 6eln 



2e(3-5e) 4(2-3e) ^ 2(l-e) 



1-e 



1-C (1-0^ 



'l-e)H''{ 6eln 



2(3-5e) ^ 4(2-3e) 2(l-e) 



l-e ■ 1 + C (l + O^J'^^^^^ 

These expressions are used in section [7] to calculate the trace of the stress energy tensor from 
the coindident propagator. 



Appendix B: Integrating the stress energy conservation 
equation 

The covariant conservation of the stress energy tensor for quantum fluctuations in FLRW 
space-times acquires in D = 4 the form (l56l) : 

1 , , 

— Pq + Apq = pg- 3pg = -Tg . (134) 

When the tress energy trace Tg is a function of the Hubble parameter H only, we have 
d/dt = —eHd/dH, such that equation (11341) can be integrated, 

45 



If, on the other hand, Tg = H^Tg{() with pg = ifVg(^), Eq. (11341) can be recast as. 




(136) 



This can be easily integrated to yield, 

P. = H^Y^^ f ^T,(C) , P, = \{p, + Tg) . (137) 
The simple useful examples that make use of Eq. (11361) are 
Tg = toH' Pi = t^H'{-^^r^) 

Tg = t,HHn{H') P. = tiH\-^^^)[HH') + ^^^). (138) 
Useful examples which make use of the integral (11371) are, 

Tg = t,H\' p^ = -t,^^C'HC) 

Tg = UH^'HO ^ p, = UH'[-^^^)[\n{C) + \). (139) 
Here t„ (n = 0, 1, 2, 3, 4) are constants, which in general depend on e and ^ but not on time. 

The special cases when u = 3/2 and z/ = 5/2 

The expressions ( 183] - [85|) and (I106H107I) are singular in the limit when u = 3/2, indicating 
that the z/ = 3/2 case requires a special attention. We shall first consider the accelerating 
case, and then the decelerating case. Furthermore, Eqs. (I114H1151) are singular in the limit 
when e = 3/2. This singular behavour can be traced back to the u = 5/2 divergence in 
Eqs. (1981) (see the n = 1 member of the sum (llOOl) ). This special case is also considered 
below. 
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9.0.4 The special case when e < 1, z/ = 3/2 



In order to get the correct stress energy tensor for z/ = 3/2 we need to take the z/ = 3/2 hmit 
of Eqs. flTTj) and fl75|) . This amounts to adding the corresponding contributions from fl78] - [80l) 
and (|83H85|) and the (7r/2) tan(7rz/) from ([H]) and ([75]). The resuh is, 



T 



C-^ (l-60^(l-e)(2-e)g^ 
' 167r2 

(l-60^(2-e)g^ 
' 647r2 

(1-60^(2-6)//^ 
' 647r2 



6e In 



6e 

Hq a 
Ha 



In 



+ 



1 4 



Hd / 4(l-e)J 
(3-56)| 



(3-5e) 



6eln 



(3-5e) 



2e' 



(140) 
(141) 
(142) 



where Hq is defined in (1ST]) . Taking a ratio of (I14ip and (11421) gives the equation of state 
parameter for this case, 



4 2e2 
-1 + -e . 

3 6e\n[Hoa/{Ha)]- {3-5e) 

4 2 
-l + -e = ^, + -e. 



(143) 



Notice that this result agrees (in the limit when a —>■ oo) with both Wg in Eq. (182]) and with 
u 3/2 limit - or equivalently the ^ —>■ e(3 — 2e)/[6(2 — e)] limit - of relation (186]) . Hence, 
when it comes to the equation of state parameter Wg at late times, there is nothing special 
about the point u = 3/2: the curves shown in figures [3H5] are continuous at z/ = 3/2. The only 
special point is the logarithmic form of the stress energy tensor, as can be seen in Eqs. ( I140f - 
11421) . One can check that the results identical to (I140H143"]) can be obtained by calculating 
the stress energy tensor from the u = 3/2 case of the half-integer coincident propagator (1^2]) . 
representing a check of Eqs. (I140H142"]) . as well as a check of our procedure based on analytic 
continuation. Notice that the logarithms drop out in the de Sitter limit when e ^ 0, which 
is a well known one loop result. The logarithms are, however, expected to re- appear at two 
or higher loop order also in de Sitter space both in massless scalar theories [29] as well as in 
quantum gravity [TO] . 



9.0.5 The special case when e> 1, z/ = 3/2 

Just as in the accelerating case, when matching onto deceleration the limit u —>■ 3/2 appears 
singular, as can be seen from Eqs. (I106H107]) . The full expressions ([96]) and ([98]) are of 
course regular (thanks for the tan(7rz^) terms). The get the correct late time limit in this 
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case, it suffices to add the tan(7rz/) terms from (1^ and (^51) to Eqs. fll06H107]) ) and take 
the z/ — > 3/2 hmit. The results are finite and - just in the accelerating case 19.0.41 - they 
acquire logarithms: 

T, (lz^myPz£)!£!c^ (..3/2) ,144) 



(lZ^M£!cMn,C) (145) 

The equation of state parameter Wq follows trivially from these relations, 

C^oo 1 e-1 1 

This relativistic fiuid scaling is in accordance with the z/ = 3/2 limit of Eq. (llOSp . Curiously 
Wq in (11471) does not depend on e, implying that as long as the relation ^ = (3e — 2)/[6(2 — e)] 
holds, Wq = 1/3 will not depend on e. Just as above, it is worth commenting that results 
identical to (I144p can be obtained by calculating the stress energy tensor from the half- 
integer coincident propagator (H2l) by taking u = 3/2 and assuming e > 1. In the minimally 
coupled scalar case when ^ = 0, the z/ = 3/2 case corresponds to matter era, e = 3/2. In 
this case the logarithmic one-loop structure exhibited in Eqs. (I144H1461) is also well known 
in literature. 

9.0.6 The special case when u = 5/2, e = 3/2 

The last special case that requires attention is the z/ — > 5/2 limit on the boundary of stability, 
where ^ = -(e - l)/[3(2 - e)] as in Eq. ffTTT]) and e = 3/2, such that Eqs. ffimim appear 
singular. This limit corresponds to the z/ = 3/2, n = 1 term in Eq. (llOOp . Similarly as above, 
in this limit one gets a finite result when the term oc C^tt tan(7rz/) from Eq. (!98|) is included. 
In this case one reproduces the z/ = 3/2 case discussed above in Eqs. (11440 with e — > 3/2: 



-^^^I^C^ (z. = 5/2, e = 3/2) (148) 



p. — i287r^ ^ ^^^^^ 
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The equation of state parameter is then, 



1 

3 



3 21n(C) 

We have thus shown that a finite answer is obtained for all values of u. 



(151) 



Appendix C: Explicitly calculating the H^(^ contribu- 
tion 

In this appendix we shall show that the contribution to the stress-energy tensor proportional 
to H^(^ is actually ultraviolet divergent. The finite contribution proportional to if^C^ will 
therefore - after renormalization - always be undetermined, until it is fixed by a measure- 
ment. Thus, not only that our procedure to calculate Pg and Pg from the trace Tg in section 
d?]) does not determine this constant uniquely, but also that in the given model the sudden 
matching at t = t generates an infinite amount of conformal fiuctuations and therefore it 
cannot be determined uniquely. In order to show this, we use (l52l) to obtain 



{n\Too\n) + 



1 



1 



D-2 

where T = T^^. Using 0541) we then find 

1-4^ 



D-2 



]iA{x;x) 



(152) 



(0|Too|n) 



D-2-A^{D-r 



{n\T\n) = dtdi + aRoo-'^tdt) ]tA{x;x] 



(153) 



Now this equation depends also on Too and therefore has a well defined contribution propor- 
tional to H'^C^, which we could not determine before, when we only calculated the contribu- 
tion proportional to T. We shall first consider the first term of the RHS 



dtdtiA(x; x] 



1 



2^-%Vr(^; 



dkk'^-^\dt^{t,k)\' 



(154) 



Now for the purpose of this section, we shall only consider the ultraviolet divergent contri- 
bution to f ll54p . We saw in section (jl]) that away from the matching point, only those terms 
that are polynomial in k contribute to the UV divergence. Thus, using ( fT9l) . we see that we 
can write 

\dt^Pit,k)\^ UV-div X\a\' + \^\')\dMty-^uit,k))\' (155) 
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Using similar techniques that led to Eq. (152]) we find in this case 



x=x,UV ^ I ' 



dkk 



D~l 



r(z/ - i + n)T{u + i + m)T{u + i + g)r(z/ - i 



"=o 1 + 1 - ^)r(// + 1 - m)r(z/ + i - g)r(i/ + | - r) 
x{^(2-.(^-,)(l-4 + ^,(l-.)^^) 
X (^4(1 - - 4z/2)(l - 4r2 - 4//^) 

clH 

+ 8i{n -r){D-l-e) ((1 - 2n)(l - 2r) - 4(1 - 2n - 2r)u^) — 

rC 

+ {D-l-ef{{l- 2nf - Au') ((1 - 2rf - Au') ^) | 



2ik 



(1 - e)p 



-m—q 



2ik 



-n—r \ 



(1 — e)aH J n\m\q\r\ 



(156) 



where as before p = a{i)H(i). To obtain the UV divergent terms, we sum the first four 
terms of (I156p to obtain 



dtdfiA^x; x) 



-D 



x=x , UV 



2^-%'^r( 



dkk 



D-l 



11/ , 



1^2 



A;2 



- + — 9-2e + e' + 2D{D - 4) - 4(1 - e)'z/ 
/ lo V 

+ s(^-'V-i)K-)^(''"-3K' 

+ (73 + - 16L>(2 - e) - 82e + 25e^ - 4(1 - efu^) 



(157) 



where C = The powerlaw divergences {k and A; ) are automatically subtracted in 
dimensional regularization and we are thus left with the logarithmic divergence. We evaluate 
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this divergence as in section 14.11 to obtain for the 1/{D — 4) part 



(2-e)(l-60 
327r2(L) -4) 



3(e(l - 2e) - 2^(2 - e)) - (2 - e)(l - Q^C (158 



where we used the expression for u from f[T5]) . The 1/{D — A) part from the other two terms 
from the RHS of fll53p are easily evaluated, using (136|) and 

Roo = -3(1 - e)H^ 
VtdtH' = 6e'H\ ^^^^^ 

We can put all terms together and obtain for the 1/{D — 4) contribution to the RHS of (11531) 
(2-e)(l-60 



327r2(D - 4) 



3e(l - 2e - 2^(1 - 4e)) - (2 - e)(l - 60C']h\ (160) 



Thus we indeed find that the contribution comes in two parts: one proportional to H^, 
and one proportional to H^C^. As a check, we can calculate the LHS of (11531) . using the 
ultraviolet contributions we obtained for the trace in (IMI) . using the conservation equation 
we can then find, using (11351) also the divergent contribution to the energy density, since if 
(n|T|f2) oc H'^ we have, apart from a possible if^^^ contribution, that 

(niTooin) = -^^^(n|T|l^). (161) 

Using (El) we then thus find for the 1/{D-A) part of the LHS offTHHl 

(2-e)(l-60 



327r2(D - 4) 



3e(^l -2e-2^(l-4e)Hi/^ (162) 



In other words, the calculation leading to (I160p is consistent with the calculation in the rest 
of this paper, apart from the H^C,'^ term. From (11601) we thus see that there is a divergence 
oc H'^C^ I (-D— 4) contributing to the stress-energy tensor. This divergence could be subtracted 
by a counterterm of the form 

a I dP Xy/^Pri (163) 



where a is a constant and pr is the pressure of some radiation fluid, obeying = |pr- This 
fluid could for example be a photon fluid, or a scalar field fluid in thermal equilibrium. This 
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renormalization is not completely satisfactory, since it requires the addition of a new field, 
not present in the original model. However, one needs to keep in mind here that our model 
is incomplete anyway. The sudden matching is put in by hand, where in a realistic model, it 
should arise from the dynamics of fields. Moreover, in a realistic model, the matching is never 
instantaneous, which should remove this UV divergence anyway. Therefore we do not feel 
that these pathologies, arising from the sudden matching are problematic. However, given 
the fact that the present model needs a counterterm like (11631) . the undetermined finite part 
contributing to the counterterm makes any H'^('^ contribution to the stress energy tensor 
arbitrary. 
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